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This paper proposes a parabolic dish/AMTEC solar thermal power system and evaluates its overall ther¬ 
mal-electric conversion performance. The system is a combined system in which a parabolic dish solar 
collector is cascaded with an alkali metal thermal to electric converter (AMTEC) through a coupling heat 
exchanger. A separate type heat-pipe receiver is selected to isothermally transfer the solar energy from 
the collector to the AMTEC. To assess the system's overall thermal-electric conversion performance, a 
theoretical analysis has been undertaken in conjunction with a parametric investigation by varying rel¬ 
evant parameters, i.e., the average operating temperature and performance parameters associate with the 
dish collector and the AMTEC. Results show that the overall conversion efficiency of parabolic dish/ 
AMTEC system could reach up to 20.6% with a power output of 18.54 kW corresponding to an operating 
temperature of 1280 K. Moreover, it is found that the optimal condenser temperature, corresponding to 
the maximum overall efficiency, is around 600 K. This study indicates that the parabolic dish/AMTEC 
solar power system exhibits a great potential and competitiveness over other solar dish/engine systems, 
and the proposed system is a viable solar thermal power system. 

© 2009 Elsevier Ltd. All rights reserved. 


1. Introduction 

Due to energy and environment problems, as well as increasing 
demand for electricity, the conversion of solar energy into electric¬ 
ity is receiving more and more attention in recent years. Sunlight is 
the world’s largest energy source and the amount that can be read¬ 
ily accessed with existing technology greatly exceeds the world’s 
primary energy consumption. Furthermore, sunlight is free, clean, 
renewable and technically exploitable in most part of the inhabited 
earth [1 ]. Today, two major solar energy technologies, photovoltaic 
and solar thermal power (electric), are being actively developed for 
electricity production. Photovoltaic, also referred to as solar voltaic, 
is the direct conversion of solar radiation to electricity using pho¬ 
tovoltaic materials. Solar thermal power (electric), also referred to 
as concentrating solar power, is electricity production from heat 
generated by solar radiation, e.g. via heat exchangers and steam 
turbines. Restricted by the capital cost of solar panels and other is¬ 
sues, the photovoltaic technology is being increasingly challenged 
by solar thermal power technology [2,3], In recent years, a number 
of solar thermal power plants have been installed in countries such 
as the US, Europe, India, and China [4,5], 

Solar thermal power systems utilize the heat generated by con¬ 
centrating and absorbing the sun’s energy to drive a heat engine/ 
generator and produce electric power. Three generic solar thermal 
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power systems, trough, power tower, and dish/engine systems, are 
being employed for power production [6,7]. Trough systems use 
linear parabolic concentrators to focus sunlight along the focal 
lines of the collectors. In a power tower system, a field of two-axis 
tracking mirrors, called heliostats, reflects the solar energy onto a 
receiver that is mounted on top of a centrally-located tower. 
Dish/engine systems, the third type of solar thermal power system, 
comprise a parabolic dish concentrator, a thermal receiver, and a 
heat engine/generator located at the focus of the dish to generate 
power. Of the three solar thermal power technologies, trough- 
electric systems are the most mature. They have a sun concentra¬ 
tion ratio of about 75 and operate at temperatures of around 400 °C 
at an annual efficiency of about 10% [8], Power towers operate at a 
concentration ratio of around 800 and temperatures of about 
560 °C and have annual efficiencies of about 15%. Dish/engine sys¬ 
tems are characterized by high efficiency, modularity, autonomous 
operation, and an inherent hybrid capability (the ability to operate 
on either solar energy or a fossil fuel, or both). A number of ther¬ 
modynamic cycles and working fluids have been considered for 
dish/engine systems. These include Rankine cycles, using water 
or an organic working fluid; Brayton cycles, both open and closed 
cycles; and Stirling cycles. Among these, Stirling engines have been 
mainly used and developed in the US and Europe [9,10]. The dish/ 
Stirling systems seem potential to become one of the least expen¬ 
sive sources of renewable energy due to their high solar-to-electric 
conversion efficiency (29.4% reported in 1984 [11], and a new 
record of 31.25% in 2007 [12]). 
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Nomenclature 



A a 

aperture area of dish concentrator, m 2 

0, 

useful energy collected, W 

Ac 

entrance aperture area of receiver, m 2 

Q rad 

radiation loss, W/m 2 

Ae 

electrode area on one BASE element, m 2 

Rb 

BASE ionic resistance, LI m 2 

A w 

cavity internal area of receiver, m 2 


the universal gas constant, 8.314J/(mol K) 

B 

temperature-independent charge-exchange coefficient, 

^int 

specific internal resistance per BASE element, Q m 2 


AK 1/2 Pa- 1 m- 2 

T a 

ambient temperature, K 

C 

geometric concentration ratio 

Tb 

temperature of BASE, K 

C pl 

heat capacity of liquid sodium, J/(kg K) 

Ted 

condenser temperature, K 

D 

diameter of the cavity, m 

T ev 

the evaporative temperature in the coupling heat ex- 

F 

Faraday’s constant, 96,485 C/mol 


changer, K 

Ib 

coefficient, F/(R g T B ), V -1 

Tev.hp 

evaporative temperature in the heat-pipe absorber, I< 

G 

dimensionless geometric pressure loss coefficient 

ATftp_ ev 

difference between T evhp and T ev , I< 

Cr L 

Grashof number based on length L 

T w 

average operating wall temperature in the cavity, I< 

g 

gravitational acceleration, 9.807 m/s 2 

A T w -h P 

difference between T w and T evhp , K 

H 

height of the cavity, m 

v° c 

open-circuit electrochemical potential per BASE ele- 


convective heat transfer coefficient, W/(m 2 K) 


ment, V 

hfg 

latent heat of vaporization of sodium, J/kg 

V cc 

v O 

effective electromotive force per BASE element, V 

I 

current through each BASE element, A 

z 

geometric factor 

Is 

direct normal insolation per unit of collector area, W/m 2 



J 

electrode current density, J = IIA E , A/m 2 

Greek symbols 

Jf 

charge-exchange current densities, A/m 2 

P 

volumetric expansion coefficient of the ambient air, I< 1 

Jsat 

saturated equilibrium charge-exchange current density, 

X 

factor of un-shading 


A/m 2 

V 

kinematic viscosity of the ambient air, m 2 /s 

k 

thermal conductivity of the ambient air, W/m I< 

£b 

radiative emissivity of BASE 

L 

characteristic dimension of cavity, m 

e c 

cavity surface emittance 

M 

molecular weight of sodium, 23 g/mol 

£cd 

radiative emissivity of condenser surface 

Nu l 

Nusselt number based on length L 

£eff 

effective infrared emittance of cavity 

Pa 

pressure of sodium vapor at the BASE/anode interface, 

(a 

concentration and charge-exchange polarization over- 


Pa 


potential at the anode, V 

poc 

1 a 

sodium vapor pressure at the BASE/anode interface in 

£e 

concentration and charge-exchange polarization over- 


open-circuit condition, Pa 


potential at the cathode, V 

K c 

pressure of sodium vapor at the BASE/cathode interface, 

0 

angle of incidence, degree 


Pa 

t/AMTEC 

thermal-electric conversion efficiency of the AMTEC 

p° c c 

sodium vapor pressure at the BASE/cathode interface in 

>le 

solar collector thermal efficiency 


open-circuit condition, Pa 

Pc 

optical efficiency 

A Pei 

total vapor pressure losses on the cathode side of the 

^overall 

overall conversion efficiency of the solar thermal power 


AMTEC, Pa 


system 

Tsat(Tev) 

saturation vapor pressure at the evaporator tempera- 

'7r 

receiver thermal efficiency 


ture, Pa 

P 

dish reflectance 

Ql 

heat losses from the receiver to the surroundings, W 

PB 

ionic resistivity of the BASE, Q. m 

Qfc 

convective heat loss through the receiver aperture, W 

TOC 

transmittance-absorptance product 

Qlk 

conductive heat loss from receiver, W 

y 

intercept factor of receiver 

Qio 

optical loss from the collector, W 

S B 

thickness of the BASE element, m 

Qlr 

radiative heat loss through the receiver aperture, W 

G 

Stefan-Boltzmann constant, 5.67 x lCT 8 W/(m 2 K 4 ) 

Or 

radiant solar energy falling on the receiver, W 

V 

tilt angle of cavity, radian 

a 

solar energy incident on the dish concentrator aperture, 

W 




However, dish/Stirling systems are quite heavy with very high 
costs. Also, because of the high pressure (>20 MPa) and high tem¬ 
perature (>700 °C), the engine are expensive to make [13]. In addi¬ 
tion, the engine in a dish/engine system converts heat to 
mechanical power in a manner similar to conventional engines; 
compressing a working fluid when it is cold, heating the com¬ 
pressed working fluid, and then expanding it through a turbine 
or with a piston to produce work. Therefore, the motion of machin¬ 
ery parts is inevitable during the operation of the dish/engine sys¬ 
tem. Moving parts are also invariably associated with problems of 
wear and tear due to the friction. Although these can be reduced 
with proper lubrication, there is virtually no way they can be elim¬ 
inated. Moreover, lubrication comes with its attendant problems 
like oil leakages, sealant problems, and possibly undesirable chem¬ 
ical interaction with other components. In addition, moving parts 
and lubricant movement (if liquid) will give rise to problems re¬ 


lated to the dynamic stability of the structure. Noise and vibration 
may also be a problem. 

A way to alleviate the problems discussed above is using a new 
power conversion unit (PCU)—alkali metal thermal to electric con¬ 
verter (AMTEC) with dishes. The AMTEC is a thermally regenera¬ 
tive electrochemical device for the direct conversion of heat to 
electrical energy [14]. It has many inherent advantages over con¬ 
ventional forms of electric generation, such as high efficiency 
(20-40%) [15], high density (up to 0.5 kW/kg) [16], absence of 
moving parts, good reliability, maintenance-free, modular design, 
competitive manufacturing costs, and compatible with many heat 
and fuel sources, including radioisotope, concentrated solar, exter¬ 
nal combustion, and nuclear reactor etc. Especially, the operating 
temperature (900-1300 K) and rated output power (several to 
50 kW) of the AMTEC perfectly matches the parabolic dish collec¬ 
tor system. 
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Furthermore, from the viewpoint of economic performance, 
AMTEC might provide significant cost savings for solar thermal 
to electricity generation. Take the dish/Stirling system as an exam¬ 
ple for comparison; the costs for a Stirling Energy Systems’ (SES) 
dish/Stirling system with rated output of 25 kW are high at about 
10,000 $/kW [9], And the PCU in the SES' dish/Stirling system costs 
33%, i.e., 3300 $/l<W. As reported in [17], the cost of AMTEC was 
estimated to be in the range between 1000 $/l<W and 2000 $/kW, 
which is 2300-1300 $/l<W lower than that of PCU in the 
SES’ dish/Stirling system. This preliminary economic comparison 
shows that adopting AMTEC as a PCU is expected to create more 
profit. 

Therefore, a parabolic dish/AMTEC solar thermal power system 
comprised of a parabolic dish solar collector and an AMTEC is put 
forward in this paper. These two subsystems are thermally cou¬ 
pled by a heat exchanger, where the useful energy collected in 
the solar collector subsystem is delivered to the AMTEC subsys¬ 
tem. The parabolic dish solar collector system consists of a dish 
concentrator, and a receiver located at the focal point of the dish. 
Existing literatures have shown that both the measured and calcu¬ 
lated total heat losses from the receiver, including conductive, con¬ 
vective, and radiative losses at higher operating temperatures are 
much more than those at lower operating temperatures [18,19]. 
That is to say, for a given type of collector, the higher the operating 
temperature, the lower thermal efficiency of the collector. On the 
other hand, for the AMTEC, the higher the temperature at which 
heat is supplied, the higher the conversion efficiency of the power 
cycle [20]. Hence, the overall conversion efficiency of the parabolic 
dish/AMTEC solar thermal power system is related to the operat¬ 
ing temperature of the coupling heat exchanger, and an optimal 
operating temperature of the coupling heat exchanger should 
exist. 

In this paper, a study on the thermal-electric conversion perfor¬ 
mance of the proposed system is presented. And this paper focuses 
on investigating the optimal operating temperature corresponding 
to the peak efficiency of the dish/AMTEC system, in conjunction 
with a parametric analysis on the overall conversion efficiency. 
This analysis is developed in four successive steps: first a simplified 
system is demonstrated, then a theoretical analysis model for cal¬ 
culating the system overall conversion efficiency is developed, 
afterwards validation of the analysis model is performed, and final¬ 
ly a parameter investigation is conducted. 


2. The description of the parabolic dish/AMTEC solar thermal 
power system 

As indicated in Fig. 1, the dish/AMTEC solar thermal power sys¬ 
tem comprises a parabolic dish solar collector subsystem and a 
power conversion subsystem—AMTEC. The two subsystems are 
thermally coupled by a heat exchanger, where the useful energy 
collected in the solar collector subsystem is delivered to the AM¬ 
TEC subsystem. 

The parabolic dish solar collector system is a point-focus collec¬ 
tor that tracks the sun in two axes. The size of the dish concentra¬ 
tor is determined by the AMTEC output power requirement. A 
separate type high temperature heat-pipe receiver is employed in 
this proposed system. Heat-pipe receiver technology has demon¬ 
strated significant performance enhancements to an already effi¬ 
cient dish/Stirling power conversion module [21]. In the 
meantime, the orientation of absorber (evaporator) of this separate 
type heat-pipe receiver is adjustable to efficiently receive solar 
rays from the dish during its rotation, while the condenser can 
be kept static. This characteristic of the receiver makes the config¬ 
uration more flexible. As shown in Fig. 1, the receiver has three 
components, including a cylinder cavity, an absorber (evaporator), 


and a condenser which also operates as the coupling heat exchan¬ 
ger. It is obvious that the coupling heat exchanger operates (func¬ 
tions) as both the condenser of the separate type heat-pipe 
receiver and the evaporator of the AMTEC. The diameter of the cyl¬ 
inder cavity is larger than its length so as to improve the absorp¬ 
tion of solar energy and reduce parasitic heat losses. Sodium 
metal is employed as working fluid to isothermally transfer heat 
from the heat-pipe receiver to the AMTEC through the coupling 
heat exchanger. This results in a uniform temperature in the cou¬ 
pling heat exchanger, thereby enabling a higher working tempera¬ 
ture for a given material and therefore higher AMTEC efficiency. In 
order to achieve the temperatures required (900-1300 I<) to effi¬ 
ciently convert thermal energy into electricity, the concentration 
ratio of the collector is typically more than 2000. To satisfy such 
a high operating temperature, structural materials should be care¬ 
fully selected. Experimental evidence has suggested that up to 
800-850 °C, stainless steel is an adequate material, while for high¬ 
er temperatures, up to 1200 °C, refractory metals (niobium, molyb¬ 
denum, tantalum), or refractory alloys, such as niobium alloys Nb- 
lZr(l%Zr) and molybdenum-rhenium (Mo-41%Re) should be used 
[17,22], 

The AMTEC subsystem consists of numbers of p"-alumina solid 
electrolyte (BASE) elements connected to improve the power out¬ 
put, because one BASE element provides relatively low voltage 
(~0.7 V). The BASE elements divide the AMTEC into two regions, 
which are the high temperature region with an evaporator and 
the low temperature region with a condenser. There are two types 
of evaporator applied in the AMTEC system. One is tube evapora¬ 
tor, which may employ an electromagnetic pump to circulate the 
working fluid from the condenser to the evaporator, as shown in 
Fig. la. The other type is shown in Fig. lb, a capillary evaporator 
based on the principle of the capillary pumped loop (CPL) is 
employed. 

The operation principle of the system is as follows: The dish¬ 
shaped concentrator reflects and concentrates incoming direct 
normal insolation to a receiver located at the focal point of the 
dish. Inside the receiver, the liquid sodium metal is vaporized 
by the intense heat created by the concentrated sunlight on the 
back surface of the cavity, on which a distributed capillary wick 
is covered. The generated sodium vapor flows up and condenses 
in the coupling heat exchanger, where it gives out the evapora¬ 
tion latent heat to the working fluid (sodium or potassium, here 
sodium is adopted as an example) in the AMTEC. The resulting li¬ 
quid sodium of the heat-pipe receiver is returned to the absorber 
(evaporator) by electromagnetic pump or passively by gravity, 
which can be achieved by the height difference between the cou¬ 
pling heat exchanger and receiver, and capillary forces in a dis¬ 
tributed wick. 

On the other hand, the liquid sodium in the AMTEC subsystem 
absorbs the heat released by the sodium vapor in the heat-pipe 
receiver loop, and becomes high temperature and high pressure 
vapor, and then is introduced into the AMTEC’s hot region. Due 
to the thermodynamic potential across the BASE, ionization of so¬ 
dium metal occurs at the anode/BASE interface. The sodium ions 
diffuse through the BASE to the cathode due to the pressure dif¬ 
ference across the BASE. The electrons circulate through the 
external load producing electrical work and then reach the cath- 
ode/BASE interface where they recombine with the sodium ions. 
The neutral sodium then reaches the condenser surface, where 
it gives up the latent heat of condensation to become liquid. 
The resulting liquid metal is circulated back by an electromag¬ 
netic pump (Fig. la) or capillary action (Fig. lb) to the coupling 
heat exchanger (the evaporator of the AMTEC). The cycle is thus 
repeated, and the AMTEC converts the work of isothermal expan¬ 
sion of sodium vapor in BASE elements directly to electric power 
[16]. 
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parabola 



A-A 


wall 


Fig. 1 . The parabolic dish/AMTEC solar thermal power system; (a) the combined system based on electromagnetic pump; (b) the combined system based on capillary 
pumped loop. 


3. Thermal-electric conversion efficiency of the parabolic dish/ 
AMTEC solar thermal power system 

The overall conversion efficiency of the solar thermal power 
system is defined as the ratio of the output electric power of the 
AMTEC to the energy incident on the dish concentrator aperture 
of the solar collector. It is the product of the solar collector thermal 
efficiency )j c and the thermal-electric conversion efficiency of the 
AMTEC >7amtec. 

t] overall = '/c'/AMTEC (1) 

The thermal efficiency of the parabolic dish solar collector sub¬ 
system and the thermal-electric conversion efficiency of AMTEC 
subsystem are analyzed as follows, respectively. 


3.1. Thermal efficiency of the parabolic dish solar collector subsystem 

The collector thermal efficiency is defined as the ratio of the 
useful energy delivered to the energy incident on the concentrator 
aperture [1], Namely, 

tic = 0-u/Q.s (2) 

Consider that the dish concentrator has an aperture area A a and 
receives solar radiation at the rate Qs from the sun, as shown in 
Fig. 2. The net solar heat transferred Qj is proportional to A a , and 
the direct normal insolation (DNI) per unit of collector area l s , 
which varies with geographical position on the earth, the orienta¬ 
tion of the dish concentrator, meteorological conditions and the 
time of day. In the present analysis, l s is assumed to be constant 
and the system is in steady state, i.e. 
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Qs = 'A (3) 

Under steady state conditions, the useful heat delivered by a so¬ 
lar collector system is equal to the energy absorbed by the heat 
transfer fluid, which is determined by the radiant solar energy fall¬ 
ing on the receiver minus the direct or indirect heat losses from the 
receiver to the surroundings. That is, 

Qu = Qr - Ql (4) 

The radiation falling on the receiver Q,- is a function of the opti¬ 
cal efficiency (/„, which is defined as the ratio of the energy falling 
on the receiver to the energy incident on the concentrator’s aper¬ 
ture. And the receiver efficiency i/ r is defined as the ratio of the use¬ 
ful energy delivered to the energy falling on the receiver [1 ]. That is 

>io = Q-r/Q-s (5) 

Vr = Q-u/Q-v ( 6 ) 

Combining Eqs. (2), (4)—(6), the collector thermal efficiency can 
be written as 



Qr 

Qs 



Or - Ql 

Qr 




= lo¬ 


th 

Qs 



( 7 ) 


where Qj is the total heat loss rate of the receiver. In order to deter¬ 
mine the thermal efficiency of the collector, the optical efficiency r\ 0 
and the total heat loss rate of the receiver Q] should be obtained. 


3.1.1. The optical efficiency 

The optical efficiency >]„ depends on the optical properties of the 
materials involved (e.g. the reflectance of dish and the optical 
properties of the receiver glazing etc.), the geometry of the collec¬ 
tor, and the various imperfections arising from the construction of 
the collector. It can be analyzed by identification of the different 
losses mechanisms [23,24]. 

• Cosine loss—quotient of total reflective area and its projected 
area, as seen from the sun. 

• Shading loss—part of the reflective area of the dish is shad¬ 
owed by the receiver, which is located at the focal point of 
the dish. But this shading can be kept a minimum if the dish 


aperture size is much larger than the receiver, generally, less 
than 1% shading [23], 

• Reflectivity loss—quotient of reflected energy and energy 
impinging on the reflective dish surface. A clean mirror made 
of low-iron glass with a silver back-coat should provide a 
reflectivity of 90-94% [23]. 

• Transmission-absorption loss—energy is lost through trans¬ 
mission of the reflected sunlight in the air from dish to the 
receiver, causing about 2-4% loss [23]; on the other hand, 
fraction of reflected radiation hitting the receiver cavity sur¬ 
face is absorbed. An ideally insulated, isothermal receiver, 
without a transparent window at the receiver entrance, can 
be considered as a blackbody with the effective absorptance 
of one [25,26]. In this study, a transparent window is placed 
at the receiver entrance to block off dust and wind, hence the 
effective absorptance of the receiver is assumed to have 
about 5-8% loss. Namely, these would cause a total of 7- 
12% loss. 

• Spillage—fraction of radiation arriving outside the entrance 
aperture of receiver may cause about 1-3% additional loss [9], 

Therefore, the following equation form can be used to perform 
an approximated optical efficiency analysis, 

r\ 0 = Xpvxy cos(0) (8) 

where X is the factor of un-shading, p is dish reflectance, xa is trans- 
mittance-absorptance product, y is the intercept factor of receiver, 
which is defined as the ratio of the energy intercepted by the recei¬ 
ver to the energy reflected by the focusing device, i.e. parabola dish, 
0 is angle of incidence. As the solar parabolic dish concentrator 
maintains its optical axis always pointing directly towards the 
sun to reflect the beam, which means the incidence angle of solar 
beam into the dish is zero degree, and the cosine loss equals to zero 
[27]. The Eq. (8) can be written as 

n 0 = *P™y (9) 

Based on the analysis above, the optical efficiency used in the 
current study, including all of these losses, is then assumed as a va¬ 
lue of 0.83, which is close to the representative value of 0.85 re¬ 
ported in [23,24]. 

3.1.2. The total heat loss rate of the receiver Qj 

The total heat loss rate of the receiver Q/ includes three contri¬ 
butions: (i) conductive heat loss from receiver, Qj k , (ii) convective 
heat loss through the receiver aperture, Qj c , and (iii) radiative heat 
loss through the receiver aperture, Qj r . The total heat loss rate (fi 
can be expressed as: 

Q/= Q«f + Q/c + Q/r (10) 

The outer surface of the receiver is usually covered with a thick 
blanket of opaque insulation to reduce the conductive heat loss. 
Research has showed that the conductive loss is normally insignif¬ 
icant compared to the convection and radiation losses [18,28]. 
Therefore, in this study, the outer walls of the receiver are assumed 
to be adiabatic and the conductive heat loss from receiver Qj k = 0. 

Among the various modes of heat loss mentioned above, convec¬ 
tion is the most complicated phenomenon and yet also a major con¬ 
tributor of the total energy loss. Hence, its characteristic has been 
extensively investigated so as to find out effective measurements 
for the improvement of system efficiency. Many authors have pro¬ 
posed models to predict convection heat loss of cavity receivers, 
such as Le Quere et al. model [29], Koenig and Marvin model [30], 
Siebers and Kraabel model [31], Clausing model [32], modified 
Clausing model [18], Stine and McDonald model [33], modified 
Stine and McDonald model [18], Paitoonsurikarn et al. model [34] 
etc. It should be noted that most models were established to 
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evaluate the natural convection heat loss, while the effect of wind 
was not considered. For the models mentioned above, Clausing 
model shows the closest prediction for downward-facing cavity 
receivers despite its original use for bigger-scale central receivers 
[35]. However, Clausing model is inconvenient to use because it is 
implicit. The simpler explicit correlation, modified Stine and 
McDonald model gave slightly better results than the implicit mod¬ 
ified Clausing model [18]. Consider that the receiver in this study 
has similar geometry and operating conditions with the receiver re¬ 
searched by Stine and McDonald, the modified Stine and McDonald 
model is selected to estimate the receiver’s convection heat loss 
while the influence of wind is neglected. Because in this study, a 
transparent window is installed at the receiver aperture to block 
off dust from the air, which also prevents forced convection flow 
to some extent, thus natural convection heat loss becomes the dom¬ 
inant loss. The modified relation of the angle dependent function 


h(cp) in [34] is adopted. It follows that 

Nu l = 0.\06Crl /3 (T w /T a ) OAS (4.256A c /A w ) s h((p) (11) 

where Gr L = gf(T w - T a )L 3 /v 2 (12) 

s = 0.56- 1.01 (A c /A w y /2 (13) 

h(cp ) = 1.1677- 1.0762 sin((p 0 - 8324 ) (14) 


In our case, the characteristic length of the receiver is the diam¬ 
eter of the cavity, D. All properties are taken at film temperature, 
i.e., the average of receiver’s cavity wall temperature and ambient 
temperature. 

The convective heat transfer coefficient can be expressed as: 
h c = Nuik/L (15) 

The convective heat loss from the receiver is given as: 

Q.ic = h c A w (T w - T a ) (16) 

Because the receiver of a parabolic dish collector system typi¬ 
cally has very small entrance aperture area A c compared to the re¬ 
ceiver’s internal cavity area A w , which results in less radiation loss. 
And the entrance aperture area A c is related to the geometric con¬ 
centration ratio C, which is defined as the ratio of dish concentrator 
aperture area to the receiver’s entrance aperture area. The follow¬ 
ing equation is used to evaluate the radiation loss from the receiver 
[36], 


T w — T ev + A T w _fy) + (sThp-ev (20) 

In general, the value of AT w _ hp + AT hp _ ev is less than 50 K, which 
is much smaller compared to the working temperature T w and T ev 
for the combined system. On the other hand, it is assumed that 
BASE temperature T B is equal to the evaporative temperature T ev 
when the heat dissipation from the coupling heat exchanger to 
BASE elements is neglected. Therefore, for convenience, it follows 
from above assumptions that T B = T ev = T w . 

The conversion efficiency of the AMTEC is defined as the ratio of 
net output electric power to all the required heat input, presented 
by the electrochemical work produced by the alkali metal ions’ iso¬ 
thermal expansion through the BASE elements at T B , the total heat 
absorbed by the working fluid to change from liquid sodium at T cd 
to vapor sodium at T ev , and the total parasitic heat losses, including 
the radiation loss, Joule heating and heat conduction loss of 
conductor leads etc. While neglecting the Joule heating and heat 
conduction loss of conductor leads, the consumed power of elec¬ 
tromagnetic pump (if used), and the parasitic losses are presented 
by radiation loss A £ q rad , the conversion efficiency of the AMTEC is 
calculated as [16] 

„ = _ KVo-R,nJ) _ 

,AMTEC Jln(P a /P^ c )/f B + [c p i(T b - T cd ) + hf g (T B )\]M/F + q rad 1 J 

In order to calculate the effective electromotive force per BASE ele¬ 
ment V C 0 C , the pressure of sodium vapor at the BASE/anode interface, 
P a , and BASE/cathode interface, P c c c , should be determined. Accord¬ 
ing to the above assumption, the pressure of sodium vapor at the 
BASE/anode interface (the high pressure side) is equal to the satura¬ 
tion pressure at T ev , namely, P a = P sa t(T e v )■ The pressure at the BASE/ 
cathode interface, Pf, is given by [38] 

P? = P° c c + AP cd (22) 

where P° c is the sodium vapor pressure at the BASE/cathode inter¬ 
face in open-circuit condition, given by Psat(T cd )y/T B /T cd , where 
Psat(Tcd) is the saturation pressure at the condenser temperature 
T cd , and AP cd is the total vapor pressure losses on the cathode side 
of the AMTEC, which is determined by [38] 

or j 

AP cd = -(2?iMR g T B ) 1/2j - (23) 


Qir = A c e eff a(Tl - T 4 ) 


where 


^eff — 1 / 


1 + 



(17) The effective electromotive force per BASE element V C 0 C , is given by 
[39] 

(18) \C = \7° c -{ fl + 4 (24) 


A c = A a /C (19) 

According to the results from [37], the effective infrared emit- 
tance of cavity £ eff is almost equal to 1 when the cavity surface 
emittance e c is in the range of 0.8-1 and the geometrical parameter 
(A w /A c ) of the focal receiver is lager than 5. Thus, the assumption of 
black body for calculating radiative heat loss through the receiver 
aperture is reasonable. 

3.2. The conversion efficiency of the AMTEC subsystem 

In reality, there exists a temperature difference A T w _ hp between 
the average operating wall temperature in the cavity T w and the 
evaporative temperature in the heat-pipe absorber T evhp due to 
the fact that the solar energy is continuously transferred to sodium 
in the absorber. Similarly, the evaporative temperature in the heat- 
pipe absorber T evhp is greater than the evaporative temperature in 
the coupling heat exchanger (evaporator of AMTEC power cycle) 
T ev , whose temperature difference is A T hp _ ev . Therefore, we have 


where V° c is the open-circuit electrochemical potential per BASE 
element, given by Nernst equation [16] 


yoc = Mb i n 
r 




(25) 


where P“ c is the sodium vapor pressure at the BASE/anode interface 
in open-circuit condition and is equal to the saturation vapor pres¬ 
sure at the evaporator temperature P sa t(T e v). 

The concentration and charge-exchange polarization overpo¬ 
tentials at the anode, i) a , and at the cathode, { c , used in Eq. (24) 
are given by the Butler-Volmer equation [40,41]: 

4 = exp(-af B {i) ~|jL eX p[(l - a )/ B {i] (26) 

Ji 1 i 

where the subscript “i” in Eq. (26) stands for either “a" (anode) or 
“c” (cathode). For a symmetric barrier, the electrochemical transfer 
coefficient a = 1/2, and Eq. (26) can be inverted to give the polariza¬ 
tion overpotential explicitly in terms of the current density: 






458 


S.-y. Wu et al./Applied Energy 87 (2010) 452-462 





(27) 


The overpotential, C i: increases logarithmically with the current, 
which reduces the voltage potential across the BASE. In Eq. (27), 
Cc < 0, Co > 0, and J c = ~J a for a thin BASE membrane. The charge-ex- 
change current densities J“, a measure of the effective local conduc¬ 
tance at the BASE/electrode/sodium vapor triple phase boundary, 
are related to the saturated equilibrium charge-exchange current 
density J“ t , as [40,41]: 


rex _ re> 

Ji —Jsi 


pot 

Kat(T B )] 


-1 1/2 


(28) 


Since the sodium vapor pressure at the BASE/cathode interface 
(normally P° c c <100 Pa) is several orders of magnitude smaller than 
that at the BASE/anode interface (normally P° a c > 20 kPa), the anode 
exchange-current density is much larger than the cathode ex- 
change-current density, hence the polarization overpotential on 
the anode side £ a is negligible. Experimental investigations have 
shown that the saturated equilibrium exchange-current density, 
J“ t , is essentially a function of the BASE temperature and of the 
type of electrodes. It can be expressed as [40,41] 


J sat = B 


Psat(T B ) 

VTb 


(29) 


Substituting Eq. (29) into Eq. (28), we have 

JT = B[P° i c P sx (T B )/T B } 1/2 (30) 


The specific internal resistance per BASE element R int , which is rep¬ 
resented by the BASE ionic resistance R B , has been given by [42] 

Rb = Pb^b (31) 


where p B is determined by Steinbruck et al. [43]. 

p B = 1.62 x lCT 5 T B exp(-45.5/T B ) + 1.55 x 1(T 7 T B exp(3722/T B ) 

(32) 


The radiation heat loss from the BASE to the condenser and the wall 
surface is most significant compared to other parasitic heat loss, 
and expressed as 

<Jrad-f0l -O (33) 


where Z is the geometric factor, given as [44] 


Z 


£ B 



(34) 


where e b is radiative emissivity of BASE, typically 0.9 above 1100 K. 
And e cd is the radiative emissivity of condenser surface, which is 
coated with liquid sodium. 


4. Validation of the analysis model 

As the proposed parabolic dish/AMTEC solar thermal power sys¬ 
tem is only a conceptual design and has not been actually con¬ 
structed in reality, therefore, no literature has been reported on 
the overall thermal-electric conversion performance of this new 
system. Fortunately, a few literatures have been published on the 
thermal performance of receiver and the AMTEC system. In this 
study, in order to evaluate the thermal performance of the para¬ 
bolic dish collector system, the optical efficiency has been approx¬ 
imately estimated, and the focus has been put on the calculation of 
heat loss in the receiver. For the purpose of validating the analysis 
model discussed in Section 3, the case of an electrically heated 
model receiver with temperatures ranging from 450 °C to 650 °C 
reported by Taumoefolau et al. [35] was studied. The results of 
comparison are presented in Table 1. We can see that the devia¬ 
tions vary from -9.1% to 17.6% for the natural convection heat loss 
at T w = 450 °C, from -13.4% to 23.7% at T w = 550 °C, and from 
-23.5% to 26.2% at T w = 650 °C, respectively. The relatively large 
deviations can be explained by two reasons. As indicated in [45], 
the calculated model adopted in this study is more accurate if 
the receiver has similar size to the data of actual receiver from 
which the model is derived. However, the fact is that the receiver 
experimentally studied in [35] was much smaller than the actual 
receiver size. The other reason is that the experimental results 
were obtained from an electrically heated model receiver in a con¬ 
trolled environment, which is different from the working condition 
of actual receivers. These two factors might together lead to the 
relatively large discrepancy. Nevertheless, the calculation results 
for the receiver considered in this study using the analysis model 
ought to be acceptable, because the receiver investigated in this 
study has similar size and working conditions with those of the de¬ 
rived receiver model. 

For validating the calculation model of the AMTEC subsystem, 
the case in [14] was studied. The efficiency of AMTEC f? amtec vary¬ 
ing with the BASE temperature T b and condenser temperature T cd 
was calculated using the model in Section 3. In [14], the cathode 
electrode was made of fine grain porous WRhi. 5 , characterized by 
B = 90 AK 1/2 /Pa m 2 , G e = 10, the current density of the BASE elec¬ 
trodes was 0.29 A/cm 2 . The results for //amtec varying with T B were 
obtained when T cd was 645 K. The results for //amtec varying with 
T cd were obtained when T B was 1127 K. The comparisons are pre¬ 
sented in Tables 2 and 3. It is found that //amtec varying with T B 
in this study is about 3.6-9.5% lower than the results in [14]. Sim¬ 
ilarly, variation of i) amtec with T cd is about 7.5-10.5% lower. The 
underestimations can be explained by the fact that, in [14], on 
top of BASE elements, thermal radiation shields were installed to 
reflect thermal radiation back towards the BASE elements, which 
reduced radiation heat loss in the low pressure cavity. For simpli¬ 
fication, this has been neglected in this preliminary study. How¬ 
ever, they show same trend and close results in the whole. 
Hence, the calculated model provided here can be used with 
confidence. 


Table 1 

Results comparison of convection heat loss from the receiver between the present analysis model and [35], 


<p (°) 

T w = 450 °C 



T w = 550 °C 



T w = 650 °C 



Qp.cal (W) 

Qa.exp (W) 

Error (%) 

Q«al (W) 

Qr.cxp (W) 

Error (%) 

Qccal (W) 

Qc.exp (W) 

Error (%) 

0 

144.1 

125.0 

15.2 

175.9 

155.8 

12.9 

208.2 

165.2 

26.2 

15 

101.4 

86.2 

17.6 

123.7 

100.0 

23.7 

146.4 

126.5 

15.7 

30 

70.9 

62.6 

13.2 

86.6 

70.8 

22.3 

102.5 

92.6 

10.7 

45 

47.2 

41.2 

14.6 

57.6 

48.6 

18.5 

68.2 

56.8 

20.1 

60 

29.6 

30.0 

-1.3 

36.2 

38.2 

-5.2 

42.8 

48.6 

-11.9 

75 

18.0 

19.8 

-9.1 

22.0 

25.4 

-13.4 

26.0 

30.8 

-15.6 

90 

12.2 

12.0 

-1.7 

14.8 

15.2 

-2.6 

17.6 

23.0 

-23.5 


Note: Subscript "c, cal” means the calculation results in this study, “c, exp” means the experimental results in [35], 
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Table 2 

Results comparison of / 7 amtec varying with T B between the present analysis model and 
[14]. 


T b (I<) 

17 AMTEC (%) in this study 

AMTEC (%) * n [ 14 ] 

Error (%) 

970 

19.0 

19.7 

-3.6 

1030 

22.5 

23.9 

-5.9 

1090 

25.0 

27.1 

-7.7 

1150 

26.6 

29.4 

-9.5 


Table 3 

Results comparison of ^amtec varying with T cd between the present analysis model 
and [14]. 


T cd (K) 

i /amtec (%) in this study 

77 AMTEC (%) in [14] 

Error (%) 

500 

25.9 

28.9 

-10.4 

550 

26.2 

29.1 

- 10.0 

600 

26.3 

29.4 

-10.5 

650 

26.0 

28.9 

- 10.0 

700 

24.8 

26.8 

-7.5 


5. Results and discussions 

To demonstrate the effects of varying parameters on the overall 
thermal-electric conversion efficiency of the parabolic dish/AMTEC 
system, numerical calculations were conducted according to the 
analysis in Section 3. Note that all the results presented here are 
based on the main parameters of the combined system summa¬ 
rized in Table 4. 

The calculated effects of varying operating temperature T w on 
tjc , ly amtec, and fjoveraii of the parabolic dish/AMTEC system are 
shown in Fig. 3. As expected, for the solar collector subsystem, 
the higher the operating temperature T w , the lower the collector 
efficiency iy c , and iy c decreases almost linearly with increasing T w . 
On the other hand, ^amtec increases evidently with increasing T w , 
which is also the operating temperature of AMTEC subsystem. 
However, when T w increases to a relatively large value, the increas¬ 
ing rate of ly amtec begins to level off. As a result, these two opposite 
effects cause ); OV eraii of the parabolic dish/AMTEC system to peak at 
20.6% when the operating temperature T w = 1280 K for the given 
conditions. However, the efficiency curve near the peak is rela¬ 
tively flat, thus it is preferable to operate on a slightly lower tem¬ 
perature, for example, T W =1200K, at the expense of a small 
decrease in )/ OV eraii (only 0.5% decrease of >yoverall, from 20.6% for 
T w = 1280 K to 20.1% for T w = 1200 K). 

As shown in Figs. 4 and 5, calculated results are presented to as¬ 
sess the effects of geometric concentration ratios C and tilting an¬ 
gle of cavity <p on i/ OV eraii for the parabolic dish/AMTEC system. 
Fig. 4 shows the )y OV eraii versus the operating temperature T w for 
various C. It is obvious that the iy ove rail increases with the increasing 
of C. And the peak ^overall increases from 19.8% to 21.2% when C in¬ 
creases from 2000 to 5000. However, it is worthy to note that the 


Table 4 

Main parameters of the combined system. 


Parameter 

Symbol 

Value 

Dish concentrator aperture area 

A, 

100 m 2 

Normal direct insolation (DNI) per unit of collector 

Is 

900 W / 

area 


m 2 

Factor of un-shading 

A 

0.99 

Reflectivity of dish surface 

p 

0.94 

Transmittance-absorptance product 

TOC 

0.9 

Intercept factor of receiver 

y 

0.99 

Aspect ratio of cavity 

H/D 

0.25 

Radiative emissivity of cavity 

S c 

0.9 

Radiative emissivity of BASE 

8 b 

0.9 

Radiative emissivity of condenser surface 

£cd 

0.02 

Ambient temperature 

Ta 

298 K 




effect of C is negligible when T w is relatively low. Conversely, the 
influence of C becomes more significant when T w is higher, due 
to that the parasitic heat losses from the receiver aperture is dom¬ 
inant at a high value of T w . Fig. 5 describes f/ 0 veraii for the parabolic 
dish/AMTEC system as function of T w , when cp = 0°, 30°, 60°, 90°, 
respectively. It can be observed that the ^overall increases with 
the increasing of <y>. The i/ OV erail is the smallest when </> = 0°, while 
the ^overall reaches the greatest value when q> = 90°. This trend is 
to be expected since as the inclination of receiver increases, more 
high temperature buoyant air remains stagnant within the cavity 
of the receiver. The maximum convection loss occurs at <y> = 0° 
where it corresponds to the peak iy OV eraii of 18.9%, and the convec¬ 
tion loss reduces to a minimum at <p = 90° where the corresponding 
peak lyoverall is 20.6%. Similar to Fig. 4, the influence of <p on /? ove rail 
becomes more considerable with the increase of T w . 

Fig. 6 illustrates the variations of fy OV eraii with T w for different 
values of T cd . All the curves show similar dependence on T w for 
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7W(K) 

Fig. 6. Variations of i/ OV eraii with T„ having T cd as a parameter. 


the Coverall of the parabolic dish/AMTEC system. It is evident that 
the i/overaii increases first, and then drops down as the T cd varies 
from 400 K to 700 K. The peak values for (/overall varies, as the 
curves show that the peak t/overaii is 20.1%, 20.4%, 20.6%, and 
19.2% when T cd is 400 K, 500 K, 600 K and 700 K, respectively. 
The present results suggest that there exists an optimum value 
of T cd for the maximum (/overall- This can be seen more clearly in 
Fig. 7, where each curve in this figure has a peak point that repre¬ 
sents the maximum f/ OV eraii corresponding to an optimal value of 
T cd . This is expected, because radiation heat loss from the BASE ele¬ 
ments to the condenser walls of the AMTEC decreases as the con¬ 
denser temperature increases. Meanwhile, an increase in 
condenser temperature causes the voltage and the electrical power 
output of the AMTEC to decrease. The combined effects make the 
>7overaii for the parabolic dish/AMTEC system peak at around 
T cd = 600 K, as shown in Fig. 7. 



However, when T cd is less than 600 K, the increment tendency of 
Coverall with T ci is quite slow. This is due to the fact that, below 
T cd = 600 K, the condenser saturation pressure is very small and 
essentially negligible compared to the pressure drop in the low- 
pressure region of the AMTEC. When T cd increases, the pressure 
at the BASE/cathode interface increases slowly, thus the voltage 
and power output of the AMTEC decreases slowly, causing a 
slow-motion of f/ OV eraii increment. While T cd is greater than 600 K, 
a sharp drop in f/ OV eraii with T cd is observed, which is explained as 
follows: above 600 K, the condenser saturation pressure increases 
exponentially with T cd , and the pressure at the BASE/cathode inter¬ 
face increases rapidly with increasing T cd , causing the voltage and 
the electrical power output of the AMTEC to decrease. Therefore, 
the decrease in heat losses and heat input to the converter as T cd 
increases, in conjunction with the decrease in AMTEC electrical 
power output, leads to the different varying tendency of (/overall 
with T cd for the solar thermal power system when T cd is less or 
greater than 600 K, as shown in Fig. 7. It should be pointed out that 
the efficiency curve near the peak is very flat, thus for the purpose 
of recycling the exhaust heat from the condenser (including the 
radiation heat loss and the latent heat released by the sodium va¬ 
por), a relatively high condenser temperature is recommended. In 
addition, it is interesting to note that the optimal value of T cd is al¬ 
ways near 600 K when T w varies, indicating that the optimal value 
of T cd is not strongly related to T w . 

The influence of dimensionless geometric pressure loss coeffi¬ 
cient G on (/overall for the parabolic dish/AMTEC system is investi¬ 
gated in Fig. 8. An increase in G results in a decrease in ;/ OV eraii. 
which is contributed to the fact that increasing dimensionless geo¬ 
metric pressure loss coefficient G means more concentration loss in 
the cathode sides of the BASE elements, which causes decreasing in 
the output of electrical power. It also can be seen that the (/overall 
decreases dramatically with the increasing of G when G is of small 
values. When G is larger, the (/overall continues to decrease, but with 
slower decreasing rate. 

The results are presented in Fig. 9 for illustrating the impact of 
the temperature-independent charge-exchange coefficient B on 
i/overaii for the parabolic dish/AMTEC system. It is shown that the 
'/overall gradually increases with the increasing of B. This trend is 
not doubtful, because the temperature-independent charge-ex- 
change coefficient B is a measurement of the contact between 
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electrode and BASE. However, the effect of B becomes negligible 
when the operating temperature T w is relatively high, which is un¬ 
like the trend in Figs. 4 and 5. A larger value of B means less internal 
charge-exchange polarization losses, thus higher value of //overall- 
This can be met by high performance electrodes, such as Rh x W 
refractory electrodes (where x varies between 1 and 2) [46,47]. 

The predicted effects of varying current density J on //overall for 
the parabolic dish/AMTEC system are presented in Fig. 10. It is 
found that the variations of //overall with T w for various J are quite 
different from those figures shown above. The curves for varied J 
in this figure cross to each other. When T w is relatively 
small (<1100 K), the lower the current density J, the higher the 
//overall- On the contraries, when T w > 1100 K, the higher the current 
density J, the higher the //overall- 

Finally, Table 5 has been presented to summarize the findings of 
the parametric analysis. 



T W (K) 

Fig. 10. Variations of ^overall with T w having J as a parameter. 


6. Conclusions 

In this study, a new parabolic dish/AMTEC solar thermal power 
system is proposed. The proposed system is an attractive solar 
thermal power system that has many advantages over the existing 
dish/engine systems. Moreover, this new system has employed a 
separate type heat-pipe receiver, in which the solar energy is iso- 
thermally transferred from the dish to the AMTEC. This type of re¬ 
ceiver makes it convenient and flexible to lay out the components 
of system. The thermal-electric conversion performance of this 
combined system has been theoretically investigated by modeling 
and calculating its overall conversion efficiency. Results show that 
the overall thermal-electric conversion efficiency of parabolic 
dish/AMTEC system peaks at 20.6% with power output of 
18.54 kW corresponding to an operating temperature of 1280 K. 
Considering the passive nature of the system, the parabolic dish/ 
AMTEC system has the potential to become one of the least expen¬ 
sive sources of renewable energy. 

The relationship between the overall thermal-electric conver¬ 
sion efficiency of the parabolic dish/AMTEC system and various 
parameters has been discussed in detail. Parameters analysis 
shows that, for the solar collector subsystem, increasing both geo¬ 
metric concentration ratio and tilting angle of the receiver would 
increase the overall efficiency of the combined system. The 
overall efficiency of parabolic dish/AMTEC system can also be 
improved by increasing the temperature-independent charge- 
exchange coefficient or decreasing dimensionless geometric pres¬ 
sure loss coefficient in the AMTEC subsystem. This improvement 
can be realized by design optimizations. The varying current 
density has distinct effect on overall efficiency of the parabolic 
dish/AMTEC system, namely, a lower current density leads to sig¬ 
nificantly higher system efficiency at a relatively low operating 
temperature, while the opposite result has been obtained at a rel¬ 
atively high operating temperature. The calculated results also re¬ 
veal that an optimum value of condenser temperature exists for 
the maximum overall efficiency. Furthermore, the maximum 
overall efficiency corresponds to a relatively high condenser tem¬ 
perature (about 600 K), which means that it is practical and cost- 
effective to recycle the exhaust heat, thus to further enhance the 
system efficiency. 
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Table 5 

A summary of the parametric analysis findings. 


Parameters Influence on r/overaii 


C 

<P 

T cd 

G 

B 

J 


(o): ^overall increases and then decreases, which resulting a peak value of 20.6% at the operating temperature T w = 1280 K for the given conditions 

(+): ^overall increases with the increasing of C, and the effect is more significant when T w is higher 

(+): ^overall increases with the increasing of <p, and the effect is more considerable with the increase of T w 

(o): r\ overall increases first, and then drops down as T cd varies from 400 K to 700 K, suggesting an optimum value of T cd = 600 K for the maximum r/ OV eraii 

(-): 7/overall decreases with the increasing of G; when G is larger, the rj oveT3 \\ continues to decrease, but with slower decreasing rate 

(+): 7/overall increases with the increasing of B, but the effect of B becomes negligible when T w is relatively high 

(-): When T w < 1100 K, Overall decreases with the increasing of J 

(+): When T w > 1100 K, ^overall increases with the increasing of J 


Note: (+) means rj overall is an increasing function of the parameter. 
(-) means r/ OV eraii is an decreasing function of the parameter. 

(o) means rj ove ra ii peaks at an optimal value of the parameter. 
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